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Cluster molecular field approximations represent a substantial progress over the simple Weiss theory 
where only one spin is considered in the molecular field resulting from all the other spins. In this 
work we discuss a systematic way of improving the molecular field approximation by inserting spin 
clusters of variable sizes into a homogeneously magnetised background. The density of states of 
these spin clusters is then computed exactly. We show that the true non-classical critical exponents 
can be extracted from spin clusters treated in such a manner. For this purpose a molecular field 
finite size scaling theory is discussed and effective critical exponents are analysed. Reliable values 
of critical quantities of various Ising and Potts models are extracted from very small system sizes. 
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I. INTRODUCTION 

The molecular field approximation (MFA) [1] is the 
method of choice to gain a first, usually qualitatively cor- 
rect, understanding of phase diagrams in cases where ex- 
act results are not available. However, quantitative pre- 
dictions derived in the framework of this approximation 
compare rather poorly with results obtained by other, 
more refined, methods. For example, the molecular field 
approximation systematically overestimates the temper- 
atures at which phase transitions take place. This is of 
course due to the fact that correlations are completely 
neglected in this treatment. Furthermore, the critical ex- 
ponents, describing the behaviour of critical quantities in 
the vicinity of a critical point, usually do not agree with 
the molecular field critical exponents. Above the upper 
critical dimension, however, the molecular field predic- 
tions of critical exponents become exact. 

The shortcomings of the original MFA led Bethe [2] 
(see also Ref. [3] and [4]) to consider a more general treat- 
ment. In the Weiss theory of ferromagnetism, a single 
spin is considered in the molecular field which is caused 
by all the other spins. This is different in the Bethe ap- 
proach where a small cluster of spins, consisting of the 
central spin and its nearest neighbours, is treated ex- 
actly. The molecular field acts exclusively on the edges 
of the cluster. Thus, correlations between a certain num- 
ber of spins arc included into the calculation. This is a 
significant improvement over simple molecular field the- 
ory, yielding for example critical temperatures which are 
much closer to the correct values than those obtained 
from the Weiss theory. 

The general idea of inserting a cluster of spins into 
a homogeneously magnetised medium led to several ap- 
plications in the past. Mullcr-Krumbhaar and Binder, 
for example, proposed a Monte Carlo method where a 
molecular field acts only on the surface spins of a clus- 
ter. This field has to be determined in a self-consistent 
way during the course of the simulation [5]. In a series 
of papers, Suzuki and coworkers established a general 



theoretical framework for extracting non-classical crit- 
ical exponents from cluster-mean-field approximations 
[6-11]. They considered a systematic series of clusters 
and showed that the true critical exponents can be ob- 
tained from the amplitudes of mean-field-type singulari- 
ties. Various other generalisation schemes of the one-spin 
molecular field approximation have been discussed in the 
literature [12,13]. 

In this work we present a systematic way of improv- 
ing the molecular field approximation by considering spin 
clusters of variable sizes embedded into a homogeneously 
magnetised background. Hereby, the interactions within 
the clusters are taken into account exactly. To this end 
we compute the density of states as a function of the in- 
ternal energy and of the interaction energy with a molec- 
ular field of unit strength. This is a hard task, but when 
it has been completed everything else does no longer de- 
mand time-consuming calculations. The actually applied 
molecular field itself is determined by imposing a condi- 
tion of homogeneity in the central part of the cluster. 
This approach is used to evaluate various temperature 
dependent quantities for Ising models defined on different 
two- and three-dimensional lattices. A two-dimensional 
Ising model with nearest and next-nearest neighbour cou- 
plings has also been studied. Results obtained for three- 
state Potts models on square and triangular lattices are 
discussed, too. 

With these data we demonstrate that critical quanti- 
ties (critical temperatures as well as critical exponents) 
can be extracted reliably from the cluster molecular field 
approximation. This is already achieved by inserting very 
small clusters containing a small number of spins (the 
largest cluster studied is formed by 181 spins) into the 
molecular field. Besides analysing effective exponents we 
also show that the cluster molecular field data exhibit 
a remarkable finite size scaling behaviour which can be 
exploited for the determination of the true non-classical 
values of the critical exponents. 

The paper is organised in the following way. In the 
next Section we present our method of systematically im- 
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proving the molecular field approximation. We also dis- 
cuss how to obtain various quantities from this data. In 
Section III the molecular field finite size scaling theory is 
discussed. We thereby use both the reduced temperature 
and the reduced energy as possible scaling variables. Sec- 
tion IV presents our data. We analyse them in different 
ways and show that critical quantities can be obtained 
from small clusters in a molecular field. Finally, in the 
last Section we present our conclusions. 



II. CLUSTERS IN A MOLECULAR FIELD 

In the following we propose a systematic scheme for im- 
proving the molecular field approximation. This scheme 
combines the exact computation of the density of states 
of small clusters with a molecular field acting solely on 
the borders of the clusters. In Section IV this method 
is applied to the determination of non-classical critical 
quantities of various Ising and three-state Potts models. 
Before discussing our approach in detail we will briefly 
reiterate the reasoning leading to the molecular field self 
consistency equation for these two models. 



A. The molecular field approximation 



In zero field the Ising Hamiltonian reads 



Hr 



(1) 



where Oi = ±1 is an Ising spin on the lattice. Angular 
brackets denote sums over nearest neighbour pairs and 
Ji is the exchange constant. The Hamiltonian for one 
selected spin i, given the spin values ai+s of its z neigh- 
bours, is 



Hi = — Ji (Ji ^ <j 



i+S- 



(2) 



<5=1 



If the spin at is considered and the rest of the lattice is 
replaced by a homogeneously magnetised medium with 
magnetisation m then the spin variables <Ji + s are replaced 
by the expectation value (cri+s) = fn and Hi becomes 
Hi^MFA^i) = -zJimvi . 

At the temperature T the expectation value of <7j in 
the molecular field z Ji to is: 



tanhfz J/m/fcgT) 



(3) 



where = cxp(z Jim/k B T) , Z- = cxp(— z Jim/ k B T) 
and Z = Z + + Z- is the partition function for the two 
states of a single spin. 

The condition (<7j) = m yields the self consistency 
equation which can be solved for the temperature 



k B T{m)/Ji = 



z to 



artanh(m) 



(4) 



Eq. (4) can be inverted or solved graphically if necessary. 
From the limit to — > of (4) one finds the MFA critical 
temperature k B T c .MF aI Ji = z and for small but finite 
to one obtains 



k B T(m)/Ji = k B T cMFA jJi 



m 



1//3a 



(5) 



with the critical exponent (3m fa = 1/2. 

The derivation of the corresponding self consistency 
equation for the three-state Potts model [14] defined by 
the Hamiltonian (the spins at taking the values 1, 2, 3) 



Hp = -J P ^ S ( a i' a i) 



(6) 



with 5{ai 1 Uj) = 1 if Ui = cij and zero otherwise is some- 
how more involved, so we only quote the result [15,16] 



to 



exp(z Jprn/ksT) — 1 
cxp(zJ P m/k B T) + 2' 



(7) 



Inspection of the molecular field free energy reveals that 
at the temperature 



k B T c /J P = 



41n(2) 



(8) 



a first order transition takes place where the order param- 
eter to jumps from to the value 1 /2. This discontinuous 
behaviour in the MFA is in contrast to the fact that the 
temperature driven phase transition is continuous for the 
three-state Potts model in two dimensions. 

These highly familiar results have been repeated also 
for the purpose of analysing the main defect of the ap- 
proach. It is the lack of any correlation between the se- 
lected spin ui and all the other spins in the lattice. This 
can be partially repaired if a cluster of N spins is inserted 
into the homogeneously magnetised medium. In our ap- 
proach, as discussed in the following, interactions within 
such clusters are taken into account exactly whereas the 
effect of the magnetised background is treated in the or- 
dinary Weiss approximation. Therefore the correlations 
of the central spin with a number of neighbouring spins 
are included into the calculation thus leading to greatly 
improved results. 



B. Cluster molecular field approximation 

In this Section we present our method to calculate 
physical quantities within an extended molecular field 
approximation based on the idea of Bcthe. 

Consider a classical spin system on a lattice defined by 
the Hamiltonian 



H = ^I(a i ,a j ) 



(9) 
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where I(<Ji, aj) is the interaction energy of the two neigh- 
bouring spins (Ji and aj. For the Ising model this 
is just the negative product of the spin variables, i.e. 
Ii(<Ti,<jj) = —JiaiCFj, for the g-state Potts model (at = 
1, . . . , q) the interaction energy is given by Ip(<Ti,<Tj) = 
— Jp5(oi,(jj). Suppose one has constructed a cluster of 
size N by some method (see below) which is embedded 
in the original lattice. For this cluster the Hamiltonian 
(9) is replaced by the molecular field Hamiltonian 

H N = ^I{<ji,(Tj) +A* I(vi,°k)5{(Jk, 1) (10) 

where indicates a summation over nearest neigh- 

bour pairs within the cluster. The second term in (10) 
represents the interaction of the cluster with the homoge- 
neously magnetised background. The sum runs over all 
pairs of nearest neighbour spins where the spin a; t is from 
inside the cluster and its partner Ok is from outside. This 
is denoted by the curly brackets. For convenience the 
are set equal to 1 and their common expectation value 
(<Tfe) = \i is placed outside the summation. The Hamilto- 
nian (10) of the extended molecular field approximation 
treats the statistical mechanics of the spin system within 
the cluster exactly and the interaction of the cluster with 
the background is modelled by a molecular field that has 
to be determined in a self consistent way. In order to 
calculate physical quantities in the cluster molecular field 
approximation it is advantageous to work out the density 
of states 



Q a0t<ri {E,B) = ]T 5(H C ,E)6(H B ,B), 

rAr(<To,Ti) 

with the internal energy of the cluster 



and the interaction energy 



H B = ^ I(ai,ak)5(ak, 1) 



(11) 



(12) 



(13) 



of the cluster with the spins outside. The summation 
in (11) runs over all Tn (<to, <ji) states of the cluster with 
fixed spin value ao of the central spin and a\ of one of its z 
equivalent neighbours. For classical discrete spin models 
the density of states is just the number of configurations 
with specified values E, B, er and o\. 

To work out physical quantities within this approxi- 
mation scheme the density of states is evaluated exactly 
so that these quantities can be computed with arbitrary 
accuracy. The technical details of computing the density 
of states (11) for a given cluster size N are summarised 
in the appendix. The size of the clusters which can be 
handled exactly is limited by the available computational 
capacities. To treat bigger clusters a recently proposed, 
highly efficient Monte Carlo algorithm to compute the 



density of states could be used [17]. However, a Monte 
Carlo calculation has to be performed with great care to 
achieve the desired accuracy. 
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FIG. 1. One of the 2 ° configurations of an Ising cluster of 
25 spins. Spin up (ex* = +1) is denoted by a full circle, spin 
down (<Ji = — 1) by an open circle. Broken bonds are marked 
with a little bar. Here E — —6 and B = —10. The central 
spin ao = +1 and one of its first neighbours o\ = — 1 are 
identified by little square. This configuration contributes one 
count to fi+_(-6, -10) = 57642. 

Clusters of increasing particle numbers are selected in 
two different ways. A possible systematic series of clus- 
ters is constructed by choosing a spin to be the central 
spin and successively including shells of nearest neigh- 
bours. Hence, the first cluster of z + 1 spins is made up 
of the central spin and its z nearest neighbours. The next 
cluster consists of the first cluster and the shell of nearest 
neighbours of the border spins of the first cluster. By this 
way a new cluster is made up of the previous one together 
with the nearest neighbours of the boundary sites of the 
original cluster. This method has the advantage that as- 
semblies of the same shape are build up. Fig. 1 displays 
the cluster of 25 spins emerging from this method for a 
two-dimensional square lattice. For the topology of the 
square lattice a cluster consists of L rows of length L per- 
pendicular to the bisector of a quadrant of the plane and 
of L — 1 short rows of length L — 1 between the longer 
ones. Hence the number of spins in a cluster is given by 
N = L 2 + (L — l) 2 . The construction of successive clus- 
ters leads to a series of clusters with 5, 13, 25, 41, 61, 85, 
113, 145, and 181 spins for the two-dimensional square 
lattice which can be handled with the available compu- 
tational capacities. For the triangular lattice clusters of 
7, 19, 37, 61, 91, and 127 spins emerge. In this geometry 
the basic row and the central row have the lengths L and 
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2L — 1, respectively. However to avoid additional com- 
putational effort due to the length of the central row for 
the treatment of the triangular lattice we alternatively 
use the clusters of the square lattice by inserting addi- 
tional bonds between the sites of the rows of lengths L 
and L — 1 so that all sites have the coordination num- 
ber 6 of the triangular topology. In three dimensions 
the equivalent procedure would be restricted to only one, 
or at most two system sizes by the available storage ca- 
pacities. This imposes a less systematic approach. The 
clusters are constructed by adding successive shells of 
neighbours of the central spin. With the nearest neigh- 
bours of the central spin one obtains a cluster of 7 spins, 
with the next nearest neighbours the size is 19, then 27, 
33, 57, and so on. The cluster of 57 spins that consists 
of the central spin and includes all spins up to the 5th 
neighbours was the largest one which could be handled 
in three dimensions. This way of selecting clusters for 
the three-dimensional case has the disadvantage that the 
various assemblies do not have the same geometry. 

Once the density of states is known physical quantities 
are easily evaluated. For example the expectation value 
of the central spin is given by 

(<j ) = ]T <r Z ao , ai /Z (14) 

with the restricted partition functions 

(E,B)exp(-{E + fiB)/k B T) (15) 

E,B 

and the total partition function 

z=J2 z °°^ • ( 16 ) 

0"O,Ci 

Similarly the expectation value of the spin a\ is computed 
from (11). The unknown homogeneous magnetisation 
/i of the sites outside the cluster must be determined 
self consistently. When the expectation value (<to) of the 
central spin is calculated as a function of /i and T the 
temperature might be adjusted such that (<to) is equal to 
H-. 

(vo)(T)=n{T). (17) 

Thereby the rigid external magnetisation is moved fur- 
ther away from the central spin and an improvement of 
the plain molecular field approximation is observed [10]. 
This is the usual condition to determine the molecular 
field discussed in the literature. Nevertheless only the 
central spin is directly related to the rigid background 
magnetisation through the self-consistency equation so 
that the expectation values of the central spin and one 
of its neighbour spins are different. In our approach the 
condition 

<<7 ) (T) = (T) (18) 



is used instead. This condition, which implicitly deter- 
mines fx as a function of T, ensures that the spin in the 
centre of the cluster and its neighbouring spins have the 
same expectation value so that at least in the central 
part of the cluster the equivalence of the spins is imple- 
mented. The spins in the centre of the cluster are not 
directly coupled to the rigid background magnetisation. 
In the extended molecular field approximation based on 
condition (18), fluctuations of the expectation value of 
a spin near the border of the cluster are suppressed due 
to the direct coupling to the homogeneous background 
magnetisation whereas the central spin can develop com- 
paratively large fluctuations. The behaviour of the cen- 
tral spin in the approximation used is consequently most 
closely related to the corresponding behaviour of a spin 
in the infinite system. The same holds for the interac- 
tion energies of bonds in the centre compared to those of 
bonds near the boundary of the cluster. These aspects 
suggest that the calculation of physical quantities should 
be restricted to the central part of the cluster where the 
condition of homogeneity is rather well fulfilled. 




FIG. 2. Order parameters as function of temperature (full 
lines) as obtained from the cluster molecular field approxi- 
mation for various two-dimensional Ising models defined on a 
square lattice. System sizes from right to left: 13, 25, 41, 61, 
85, 114 145, and 181 spins. The spontaneous magnetisation 
of the infinite system is shown as a dashed line. For sys- 
tems treated in the cluster molecular field approximation the 
order parameter sets in abruptly at a pseudo-critical temper- 
ature, see Table I. At low temperatures the finite size order 
parameters follow closely the spontaneous magnetisation of 
the infinite system. 

After solving the implicit equation (18) for the molec- 
ular field all physical quantities can be evaluated using 
the density of states (11). The order parameter of the 
system is determined from (cr ). Note that in the limit 
T — > the interaction of the border spins of the cluster 
with the homogeneously magnetised background forces 
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all spin variables to be in the spin state a = 1 of the 
spins outside the cluster. For the Ising model the order 
parameter is given by 

mi (T) = (<7 ) (T) (19) 

and for the g-state Potts model one has 

mp(r) = g(g ° g )(T i ) ' 1 . (20) 

In the cluster molecular field approach the order pa- 
rameter is exactly zero at high temperatures. When 
decreasing the temperature the order parameter of the 
Ising model exhibits a sudden onset at the tempera- 
ture T c n indicating a transition to a magnetised phase. 
The abrupt appearance of a finite order parameter de- 
fines the pseudo-critical temperature of the finite sys- 
tem. Fig. 2 illustrates this temperature dependent 
behaviour of the order parameter in the case of the 
square lattice Ising model. It is instructive to compare 
the pseudo-critical temperatures resulting from equa- 
tions (17) and (18). For a cluster of TV = 145 spins 
defined on the square lattice our approach (18) leads 
to the value k B T c , N jJ I = 2.37096 (T C ^ N /T C = 1.045, 
where ksT c /Ji = 2/ln(l + \/2) is the critical tempera- 
ture of the two-dimensional Ising model defined on the 
square lattice), which should be compared to the value 
2.57524 (T C . N /T C = 1.135) obtained when using the self- 
consistency equation (17) [10]. A similar remarkable im- 
provement is also achieved for three-dimensional clus- 
ters: for N = 57 spins our pseudo-critical tempera- 
ture is fc B T c , w /J/ = 4.6627 {T CjN /T c = 1.034), much 
closer to the value 4.5115 of the infinite system than 
the value 5.0018 (T C . N /T C = 1.109) resulting from Eq. 
(17) [10]. The pseudo-critical temperatures of the Ising 
model for various lattice topologies are summarised in 
Table I and Table II for the two-dimensional and the 
three-dimensional systems, respectively. 

In Fig. 3 we show the order parameter for the three- 
state Potts model defined on a square lattice as a function 
of temperature for a series of cluster sizes. For all clusters 
a jump in the order parameter is observed, thus disclos- 
ing the discontinuous character of the phase transition 
in the MFA. This jump gets smaller when the number 
of cluster spins is increased reflecting the softening of 
the first order character of the phase transition with in- 
creasing cluster size. In the limit TV — > oo the jump 
will disappear completely in accordance with the contin- 
uous phase transition taking place in the bulk three-state 
Potts model. The temperature at which the jump occurs 
is deduced from the behaviour of the free energy. The 
values of these transition temperatures are displayed in 
Table III for the square and the triangular lattice. 

We estimate the internal energy from the bonds of the 
central spin with its z neighbour spins which yields 

u{T)= Z -Y j no^a x )Z aQ . a jZ . (21) 



TABLE I. The pseudo-critical temperatures T Cj jv of the 
two-dimensional Ising model defined on the square (sq) and 
the triangular (tr) lattice. 



N 


sq 


tr 
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145 


2.37096 


3.82574 


113 


2.38480 


3.85042 


85 


2.40294 


3.88338 


61 


2.42778 


3.92602 


41 


2.46380 


3.99458 


25 


2.52078 


4.0804 


13 


2.62534 


4.27958 


5 


2.88539 


4.72088 



TABLE II. The pseudo-critical temperatures T Ci at of the 
three-dimensional Ising model on the simple cubic (sc) lattice. 



N 


sc 


oo 


4.5115 


57 


4.66266 


33 


4.69361 


27 


4.69367 


19 


4.74872 


7 


4.93261 



TABLE III. The transition temperatures of the three-state 
Potts models defined on the square (sq) and the triangular 
(tr) lattice. 



N 


sq 


tr 


oo 


0.99497 


1.58493 


85 


1.01937 


1.62849 


61 


1.02396 


1.63619 


41 


1.03051 


1.64813 


25 


1.04074 


1.66267 


13 


1.05957 


1.69818 
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FIG. 3. Temperature dependent order parameters of 
three-state square lattice Potts models derived within the 
cluster molecular field approximation (full lines) . The exactly 
treated clusters contain (from right to left) 13, 25, 41, 61, and 
85 spins. The discontinuous character of the phase transition 
in the cluster molecular field approximation is displayed by 
the jump of the order parameter. The unstable solution for 
a system with 13 spins is shown as a dashed line. The grey 
line indicates the location of the critical point of the infinite 
system: k B T c /J P = l/ln(\/3 + l) w 0.995. Inset: free en- 
ergies obtained for a cluster with 13 spins within the cluster 
molecular field approximation. The free energy of the sta- 
ble high and low temperature phases are shown as full lines 
whereas the dashed line is the free energy of the unstable so- 
lution. The intersection of the full lines yields the transition 
temperature listed in Table III. 

From the internal energy one now gets the specific heat 

du{T) 



c{T) 



dT 



(22) 



The specific heat is also related to the entropy of the 
central spin through the relation 



c(T) = T 



ds(T) 
dT 



(23) 



This allows the computation of the free energy 

f(T) = u{T) - Ts(T) (24) 

from the knowledge of the internal energy. Solving equa- 
tion (23) for the entropy s(T) one obtains the free energy 



f(T) = u(T) -T f dr^ 
Jto t 



(25) 



with To being some arbitrarily chosen temperature rep- 
resenting the constant of integration. The free energy of 
the Potts model for the cluster with 13 spins in the vicin- 
ity of the transition point is shown in the inset of Fig. 
3. The solution of (20) that minimises the free energy is 



physically stable. This leads to the kink of the free en- 
ergy at the transition temperature again indicating the 
discontinuity of the phase transition for finite N. 

Next we consider the response of the central spin to 
an external field that acts exclusively on the central spin 
itself. In order to evaluate this zero field susceptibility 
one has to differentiate the magnetisation with respect 
to a magnetic field that acts only on the central spin 
variable ctq 1 



X(T) 



lim 

h— >o 



dm(T, h) 
dh 



(26) 



To obtain an analytic expression for the susceptibility 
of the Ising model we add a magnetic field term to the 
Hamiltonian: 



Hj 



(27) 



When evaluating the susceptibility one has to bear in 
mind that one gets a field-dependent molecular field fx = 
/x(T, h) due to this additional external field h. This has 
the consequence that (26) depends on the derivative of \x 
with respect to the magnetic field h. However, the order 
parameter m is given by the two equivalent equations 
m = ((To) and m = (<7i) so that one of them can be used 
to eliminate d^/dh from (26). Differentiating (oo) with 
respect to h one obtains the expression 



x(T,h) 



1 



<-o 2 > 



dh 



«<7 > (B) - (a B)) 



(28) 



On the other hand differentiation of (cti) yields 
1 



X(T, h) 



k B T 



^o) 2 + ^«a ) (B)-(a 1 B)) 



(29) 



Combining these two relations in order to eliminate the 
derivative of the molecular field results in the desired ex- 
pression of the zero-field susceptibility for temperatures 
below T c jy: 



X(T) 



1 



{4) W 



({a ) (B)-(<T B))((^)-(a a 1 )) 
(a B) - (o-iB) 



(30) 



The first two terms in the square brackets represent 
the ordinary fluctuation formula for the susceptibility. 
The additional third term arises from the coupling to 
the molecular field. Its denominator vanishes in the 
limit T — ► T c jq yielding a diverging susceptibility for 
T — > T c jv as expected for the molecular field approach. 
Fig. 4 displays the zero field susceptibility of the two- 
dimensional Ising model defined on the square lattice for 
different cluster sizes. 
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A. Finite size scaling for the temperature picture 

The spontaneous magnetisation of the cluster molec- 
ular field approximation displays the characteristics of 
spontaneous symmetry breaking. This allows an unam- 
biguous definition of the pseudo-critical temperature T c ,l 
of the finite system of linear extension L taken to be N 1 ^ 
where d is the dimensionality of the lattice. The spon- 
taneous magnetisation near this pseudo-critical temper- 
ature then behaves like 



k B T/J, 

FIG. 4. Zero field susceptibilities as function of tempera- 
ture obtained from the cluster molecular field treatment of 
two-dimensional square lattice Ising clusters, see Eq. (30). 
The susceptibilities diverge at the pseudo-critical tempera- 
tures. The Ising clusters contain 25, 61, 113, and 181 spins. 

For the g-state Potts model the magnetic field acts on 
the central spin only if it is in the spin state a of the sites 
outside the cluster: 



m sp , L (t') = A L (-t'f 
where the reduced temperature t' is defined as 
T - T c .l 



Hp — + Hp — hS(ao, a) 



(31) 



The susceptibility of the Potts model for temperatures 
below the pseudo-critical point of the finite cluster is then 
given by 



X(T) 



1 



[(a Q S(a Q ,cr)) - (a ) (6(<r a ,a)) 



q-1 k B T 
({ao) (B) - (aoB))({(ToS(a ,a)) 



(<t q B) - (axB) 



This expression for the Potts model has the same struc- 
ture as the corresponding expression (30) of the Ising 
model. However the susceptibility of the Potts model 
does not diverge as the molecular field displays a jump 
at the transition temperature. 



III. MOLECULAR FIELD FINITE SIZE SCALING 
THEORY 

In this Section we present a phcnomcnological finite 
size scaling theory for our cluster molecular field approx- 
imation. We consider both the reduced temperature and 
the reduced energy as scaling variables and discuss the 
relative size of the corrections to scaling. Our data are 
then analysed in the following Section using the frame- 
work of the molecular field finite size scaling to be ex- 
posed in the following subsection. 



t' = 



T 



c,L 



(33) 



(34) 



and where /3' = 0.5 is the molecular field critical expo- 
nent. For increasing clusters the influence of the molecu- 
lar field which acts only on the boundary becomes weaker 
so that the approximation eventually becomes exact in 
the limit L — > oo. In the infinite system the spontaneous 
magnetisation vanishes like 



m sp (t) = A(-t) 
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(35) 



for small reduced temperatures t which is related to the 
critical temperature T c — T c ao of the infinite system by 



t = 



The fractional shift rj, is given by [18] 



T n — T. 



(32) 



c.L 



—DL~ 



(36) 



(37) 



with the shift exponent A = 1/v where v is the critical 
exponent of the correlation length. Equation (37) and 
(34) are combined to yield 



t' = t + t l + 0{tI) 



(38) 



for large L. Equation (37) can be used to obtain an 
approximation of the transition temperature in the ther- 
modynamic limit. This is exploited in Section IV. 

For sufficiently small t' we assume that the sponta- 
neous magnetisation of the finite system of extension L 
is given by 



m sp , L (t') = LP ( f>(t'L 1 ^) 



(39) 



The function 4>{x) is universal showing no further L- 
dependence. The exponent p in (39) is obtained by taking 
the limit L — *■ oo for fixed temperature T ^ T c . In this 
limit equation (39) reduces to (35) yielding 



4>(x) 



as x 



(40) 
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The residual L dependence is cancelled out if p — — (3 /v. 
Hence the finite size scaling relation 

L^m sp , L (t') = WL 1 '") (41) 

follows. 

At this point we have to mention that an unambiguous 
definition of the linear extension of our finite clusters 
is not possible. This crucial fact was already pointed 
out in [7]. For our investigations we choose L = N 1 ^ 
but other definitions are also possible. For the three- 
dimensional systems our choice is especially problematic 
as the clusters do not have the same geometry. We will 
come back to this point in Section IV when discussing 
our numerical data. 

From equation (33) follows that the universal function 
<j> behaves like a square root for small scaling variables, 
i.e. 

^{t'L 1 '") = A+yZ-t'L 1 /" , (42) 

for t' — > and fixed L. Taking the limit t' — > and 
comparing with relation (39) the finite size scaling law 

A L = A^L^ (43) 

for the amplitudes of the magnetisation of the finite sys- 
tems is derived. 

In an analogous way we obtain a finite size scaling rela- 
tion for the susceptibility. In the molecular field approach 
the susceptibility diverges with the molecular field criti- 
cal exponent 7' = 1. We then have the scaling relation 

L-^ X L(t') = ^(t'L 1 '") , (44) 

where for small scaling variables the universal function ip 
is given by 

^(t'L 1 '") = A^p {t'L 1 '")- 1 . (45) 

One also derives, in complete analogy with Eq. (43), a 
finite size scaling law involving the finite size amplitudes 
of the susceptibility. Note that this relation was already 
used in references [7,10] to get estimates for the criti- 
cal exponent 7. All these relations show that the critical 
exponents of the infinite system are encoded in the ampli- 
tudes of the corresponding quantity of the finite system. 
It is precisely this fact which will allow the determination 
of critical exponents from the cluster molecular field ap- 
proximation. It is also worth noting that in the present 
approach the critical temperature of the infinite system 
is not needed for the determination of critical exponents. 

B. Finite size scaling for the energy picture 

The internal energy of a system is given by the expec- 
tation value (21). The expression ul = ul(T) can be 
inverted in order to obtain the temperature as a function 



of the internal energy of the system. All physical quan- 
tities can then be expressed as functions of the internal 
energy. In the following subsections this transformation 
from the temperature picture to the energy picture is dis- 
cussed for systems displaying either an algebraically or a 
logarithmically diverging specific heat in the thermody- 
namic limit. 



1. Algebraically diverging specific heat 

Consider the spontaneous magnetisation near the 
pseudo-critical point of the finite system: 

m sp , L (t') = A^L i ^ £ V^t J (l + B L t' + ...) (46) 

where a correction term to the leading behaviour with 
the L-dependent amplitude Bt, is explicitly taken into 
account. From finite size scaling theory one obtains a 
similar expression for the (intensive) zero-field specific 
heat near t' = for negative t': 

c h=0 .L{t') = aL^{l + 2b L t' + ...) (47) 

with an L-dependent amplitude of the linear correction 
term. Here it is assumed that the singularity of the spe- 
cific heat of the infinite system is characterised by a 7^ 0. 
The case of a logarithmic divergence is deferred to the 
next subsection. Note also that the specific heat in molec- 
ular field approximation does not diverge. It exhibits a 
jump singularity, hence one has the exponent a' — 0. At 
the pseudo-critical point of the finite system, i.e. t' = 0, 
equation (47) reduces to 

C L = c fc=0 ,L(0) - aL^ . (48) 

This relation can be used to estimate the true critical 
exponent a from the maxima of the specific heat of the 
various finite systems. 
As 

_ du{T) _ du dt' _ 1 du 
c h=0 .L - - wHT ~ TTlW { ' 

integration of (47) gives 

e'{t')= u L {t')-u c , L (50) 

with u Ct L being the constant of integration. This con- 
stant is just the pseudo-critical energy of the finite sys- 
tem. The expression for e'(t') can be inverted for small 
t' to give 

t> = ^L-^e' ^_L- 2 «/V 2 + . . . (51) 

J c,L« 1 c,L a 

in the asymptotic limit e' — > 0. Plugging this result in 
equation (46) one finally obtains 
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m sPtL (e') = A^L^^ 



aT, 



c,L 



the two-dimensional Ising model. In this case the ampli- 
tude Cl (i.e. the maximum) of the finite system specific 
heat varies as [20] 



(l + L- a l"-^-{B L + • • •) ( 52 ) C L = A\nL"» + £ = A' InL + £ . 



(55) 



with the renormalised exponents /3 £ = and v e = 

jz^- These renormalised critical exponents describe the 
singular behaviour of the various thermostatic quantities 
of the infinite system in the energy picture [19]. The 
reformulation also yields the finite size scaling relation 



L ' /v 'm sPiL (e') = £ ( 



(53) 



for the energy picture with a scaling variable that con- 
tains the transition temperature of the finite system. In 
a completely analogous way the reformulation can be 
done for other physical quantities like the susceptibility. 
For large systems the pseudo-critical temperature T c .l is 
close to the critical temperature of the infinite system so 
that is can be replaced by T c . However, as we are deal- 
ing with very small clusters the additional L-dependence 
in the scaling variable due to T c l can not be neglected. 
Similar to equation (37) one obtains the shift relation for 
the pseudo-critical energy: 



-D F L 



-1/Ve 



(54) 



The proposed scaling variable for the energy picture and 
the corresponding scaling relations are tested in Section 
IV. 

In equation (52) the amplitude of the linear correc- 
tion term contains the amplitudes B L and br, of the tem- 
perature picture. Moreover a prefactor L^ a / U appears 
which will lead to a suppression of the leading correction 
term compared to the corresponding expression (46) in 
the temperature picture if the critical exponent a is posi- 
tive. This suppression is determined by the prefactor to- 
gether with the L-dependence of the amplitudes Bl and 
hi,. If the amplitudes grow strongly enough with the sys- 
tem size the suppression becomes weaker for increasing 
L. Note however that a reduction of the linear correc- 
tions in the energy picture will eventually always take 
place even if this reduction becomes weaker for larger 
cluster sizes. For a system with a negative exponent the 
correction term in the energy picture is blown up relative 
to the temperature picture. The investigation of this in- 
teresting consequence in the 3d Heisenberg model as an 
example of a system with a < is left to future studies. 



2. Logarithmically diverging specific heat 

With slight modifications the above discussed trans- 
fer to the energy picture is carried out for the case of a 
logarithmically diverging specific heat c = — A\nt of the 
infinite system. A prominent example of such a system is 



O 




In L 

FIG. 5. Specific heat maxima as function of InL for small 
two-dimensional Ising clusters defined on a square lattice. 
The linear extension L is given by L = \fN where N is the 
total number of spins forming the cluster. The dashed line 
results from fitting these data to Eq. (55), with A' = 0.562 
and £ = 1.76. 

Fig. 5 shows equation (55) for the 2d Ising model in the 
extended molecular field approximation. The depicted 
linear behaviour is characterised by the slope A' = 0.562 
and the intercept £ = 1.76 with the vertical axis. For suf- 
ficiently large system sizes the L-independent constant £ 
can be neglected relative to the increasing contribution 
from the logarithmic term. As relatively small systems 
are considered in the molecular field approximation dis- 
cussed here this constant is explicitly taken into consid- 
eration. The arguments used to replace the variable t' by 
the internal energy still apply and give rise to the relation 



m S p,L(s') = A^L "e 



-£' 



{A' In L + £)T C , L 



(A' In L + 0T c .L 



(B L - b L )e' + . . . . (56) 



The linear correction term is again diminished by an ad- 
ditional L-dependent prefactor as already observed for 
systems with a > 0. This form of the magnetisation 
suggests the finite size scaling law 



L^m sp , L (e') = . 



z , L i/ Ve 



{A' In L + C)T C , L 



ClT c x 
(57) 



with a slightly modified scaling variable. This modi- 
fied scaling variable in the energy picture is investigated 
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for the two-dimensional Ising model in the next Section. 
Note however that the factor causing the suppression of 
the linear correction term of the magnetisation is now 
integrated into the scaling variable. This has the conse- 
quence that the reduction of the corrections to scaling in 
(57) are solely controlled by the relative size of the am- 
plitudes Bl and &l. Therefore it has to be expected that 
the reduction of the corrections are weaker than those in 
the case of an algebraically diverging specific heat. 



IV. CRITICAL QUANTITIES 

In this Section we extract critical quantities like crit- 
ical exponents or transition temperatures of the infinite 
system from the molecular field data of the correspond- 
ing finite systems. Apart from the molecular field finite 
size scaling relations presented in the previous Section 
the notion of the effective critical exponent provides an- 
other way to obtain critical exponents from finite system 
quantities [21]. 

The effective critical exponent f3 e ff,N of the order pa- 
rameter is defined by expressing the temperature depen- 
dent spontaneous magnetisation in the form 



m N = At p ""- N{T) 



(58) 



where t = (T c — T) /T c and T c is again the critical temper- 
ature of the infinite system. The temperature dependent 
effective exponent is obtained by differentiating relation 
(58): 



dln(mjv) 
dint 



(59) 



Fig. 6 shows the effective exponent of the two- 
dimensional Ising model for the different clusters on the 
square lattice together with the effective critical exponent 
Peff, oo °f the infinite system [22]. As the pseudo-critical 
temperatures of the finite systems are above the critical 
temperature T c of the infinite system the spontaneous 
magnetisation itin(T) is finite at T c and therefore the ef- 
fective exponent has to vanish at T c . It is also observed 
that the effective exponents of the largest clusters are 
closer to the corresponding curve of the infinite system 
than those of the small clusters. From the effective expo- 
nents estimates for the true critical exponent f3 may be 
extracted by means of linear extrapolation. To this end 
an extrapolated exponent (3 ex t,N is defined by 



[3ext,N( T )=Peff,N(T) + 



dT 



(T c - T) . (60) 



This extrapolated exponent is displayed in Fig. 7 for the 
square lattice again together with (3 ex t,oo of the infinite 
system. Note that the effective exponent depends on the 
critical temperature of the infinite lattice. For models 
with an unknown value of T r . a reliable estimate has to 



be obtained from the finite system pseudo-critical tem- 
peratures. This can be achieved by using the finite size 
scaling relation (37). 
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FIG. 6. Effective critical exponents /3 e ff,N versus temper- 
ature for different square lattice Ising clusters. The system 
sizes are the same as those in Fig. 2, the lowest line resulting 
from the cluster with 13 spins. These curves directly follow 
from the logarithmic derivative of the spontaneous magneti- 
sations displayed in Fig. 2 with respect to the temperature. 
The dashed line is the effective critical exponent of the in- 
finite system, yielding the critical exponent (3 = 1/8 in the 
limit T — > T c . 
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FIG. 7. Temperature dependent extrapolated exponents 
Pext,N(T) derived from the effective critical exponents shown 
in Fig. 6, see Eq. (60). The dashed line is again the corre- 
sponding quantity of the infinite system. The inset demon- 
strates a peculiarity of the 2d Ising model. The extrapolated 
exponent of the infinite system overshoots the asymptotic 
value 1/8 and this is also true for the cluster MFA results 
for N = 145 and 181. 

In the next Subsections we investigate our molecular 
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field data for the various systems using both finite size 
scaling relations and effective exponents. 



A. Two-dimensional Ising models 

Using the cluster molecular field approximations we 
studied two-dimensional nearest neighbour Ising models 
defined on the square and on the triangular lattices as 
well as an Ising model with both ferromagnetic nearest 
and next-nearest neighbour interactions defined on the 
square lattice. The extrapolated exponent from Eq. (60) 
allows us to obtain reliable values for the order parameter 
critical exponent (3, as shown in Fig. 7 for the square 
lattice model with only nearest neighbour interactions. 
The best estimate for f3 is obtained from the value of 
the temperature where the extrapolated exponents for 
the two largest systems (N = 145 and N = 181) still 
coincide, yielding (3 w 0.126 which should be compared 
to the exact value (3 ex — 1/8. This overestimation of 
only 0.5% is rather remarkable considering the smallness 
of the systems under consideration. Note however that 
the slight overshooting above the true value 1/8 is not 
an artifact of our approximation, it is due to the peculiar 
behaviour of the exact solution of the infinite system, see 
the inset of Fig. 7. 

A similar behaviour is also observed for the other two- 
dimensional Ising models. In these cases the overshooting 
is somewhat larger, yielding a slightly increased value for 
(3. Thus we obtain the value (3 « 0.128 for the triangular 
lattice, still in very good agreement with the exact result. 

Other critical exponents, as for example the critical 
exponent 7 of the susceptibility, might be derived in the 
same way. In the present study, however, we restrict 
ourselves to the small system sizes where the density of 
states can be computed exactly. For these systems the 
extrapolated exponents obtained from the susceptibility 
do not yet present a plateau, making an estimation of 7 
rather tedious. Here, the molecular field finite size scal- 
ing theory of Section III is an interesting alternative for 
extracting critical quantities, as discussed in the follow- 
ing. 

Let us start by probing the finite size scaling relations 
(41) and (57) of the spontaneous magnetisation in the 
temperature as well as in the energy picture. Note that 
relation (57) has to be used instead of (53) as the specific 
heat of the two-dimensional model presents a logarithmic 
divergence with a — 0. Plotting the scaled magnetisation 
as a function of the scaling variable and setting L — N 1 ^, 
d being the number of space dimensions, the best data 
collapse as derived from a least-square fit procedure or 
from a more sophisticated approach [23] yields the val- 
ues of the critical exponents. The resulting data collapses 
are shown in Fig. 8 for the nearest neighbour Ising model 
defined on the square lattice where the smallest systems 
with 5 and 13 spins have been discarded. Hereby the crit- 
ical exponents in the temperature picture take the values 



(3 = 0.127 and v = 1.01 (see Fig. 8a), while in the energy 
picture (see Fig. 8b) the best data collapse is achieved 
for (3 £ = 0.124 and v £ = 1.02. This convincingly shows 
that the finite size scaling relations derived in Section III 
are indeed fulfilled by our scaled data. Again, one has to 
stress the remarkable agreement between our estimates 
of the order parameter critical exponents and the exact 
value 1/8, based on data obtained for very small system 
sizes. 
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FIG. 8. Scaled magnetisations of Ising clusters defined on a 
square lattice, see Fig. 2, as function of different scaling vari- 
ables: (a) in the temperature picture, see Eq. (41), and (b) in 
the energy picture, see Eq. (57). The best data collapse yields 
the following values for the critical exponents: (a) (3 = 0.127, 
V = 1.01, (b) /3e = 0.124, v e = 1.02. 

Fig. 9 displays the scaled susceptibilities as function 
of the scaled reduced temperature. The values of the 
critical exponents 7 = 1.77 (to be compared with the 
exact value "f ex = 1.75) and v = 1.01 again result from 
the best collapse of the scaled data. The same can be 
done in the energy picture where we obtain j e = 1.80 
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and v e — 1.02, in reasonable agreement with the exact 
values. This shows that the finite size scaling theory 
enables us to obtain reliable values for critical exponents 
also in cases where an analysis of effective exponents does 
not succeed. 
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FIG. 9. Scaled susceptibilities obtained for square lattice 
Ising clusters in the cluster molecular field approximation, 
see Fig. 4, as function of the scaled reduced temperature. 
The values 7 = 1.77 and v — 1.01 of the critical exponents 
result from the best data collapse. 

At this stage we have to emphasise again that no 
non-universal quantities of the infinite system enter the 
molecular field finite size scaling theory. Especially, the 
exact location of the critical point is not needed. On the 
other hand, using the value of v obtained from finite size 
scaling together with the pseudo-critical temperatures of 
the three largest systems listed in Table I, a remarkable 
good estimation of the critical temperature follows from 
Eq. (37): ksTc/Ji — 2.266 in excellent agreement with 
the exact value fc B T c /J/ = 2/ln(l + y/2) w 2.269. 

Before closing this subsection, we note that almost 
identical estimates are obtained for critical quantities 
when studying the triangular model: /? = 0.127, (3 £ = 
0.125, 7 = 1.78, 7 £ = 1.80, v = 1.03, v e = 1.03, and 
UbTc/Ji = 3.631 (exact: 3.634). The discrepancies are 
somewhat larger for the square lattice clusters where the 
next-nearest neighbour interactions have also been taken 
into account. We obtain f3 = 0.120 and 7 = 1.86, for 
example, analysing systems with up to 81 spins. 



B. Three-dimensional Ising model 

When trying to extract critical quantities for the three- 
dimensional Ising model from the cluster molecular field 
approximation, we are faced with two difficulties: (1) 
only systems with a very small number of spins (< 57) 
can be handled exactly and (2) the various clusters have 
different geometries, as explained in Section II. Espe- 



cially the second point makes a finite size scaling analysis 
much more problematic, as discussed below. But let us 
start again by looking at the effective exponent derived 
from the spontaneous magnetisation. Fig. 10 shows the 
extrapolated exponent obtained for the different system 
sizes. The estimated value of 0.316 for the critical expo- 
nent /3 is slightly smaller than the expected value 0.326. 
This is due to the finite-size effects which for our small 
systems set in already at rather low temperatures. Fig. 
10 also displays an oddity due to the different cluster 
geometries: on the dashed line the extrapolated expo- 
nents of two different system sizes, N = 27 and N = 33, 
are almost indistinguishable. Looking at the raw data, 
one observes that, by coincidence, due to the difference in 
the geometries, the cluster molecular field approximation 
yields virtually identical spontaneous magnetisations for 
these two clusters. 
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FIG. 10. Extrapolated exponents (3 e xt,N versus tempera- 
ture for different three-dimensional Ising clusters defined on a 
cubic lattice. The quantities derived for systems with 27 and 
33 spins are indistinguishable (indicated by the single dashed 
line) . 

The difficulties originating from the difference in ge- 
ometry are enhanced when looking at the finite size scal- 
ing behaviour in three dimensions. In fact, the rather 
naive identification L = N 1 ^, which works very well in 
two dimensions with identical geometries, is inadequate 
when different cluster geometries are involved. Sticking 
to this definition of L, we are only able to obtain simul- 
taneously a good data collapse for all the quantities of 
interest (magnetisation and susceptibility as function of 
both reduced temperature and reduced energy) when the 
sizes TV = 27 and TV = 33 are omitted. This procedure 
then yields the following values for the critical exponents: 
f3 = 0.32, e = 0.41, 7 = 1.23, 7e = 1.45, v = 0.63 and 
v £ = 0.69. Whereas in the temperature picture the val- 
ues of the critical exponents agree nicely with the values 
given in the literature (/? = 0.326, 7 = 1.24, v = 0.63), 
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some discrepancies show up in the energy picture (for 
example, one expects 0.367 for (3 e ), even if the overall 
agreement is still very satisfactory with regard to the 
smallness of our systems. 

A possible way to circumvent some of the difficulties 
arising from the difference in cluster geometry would be 
to eliminate L altogether in Eq. (41) resp. (53) by ex- 
ploiting [7] the relation (37) resp. (54). However, this 
would introduce the necessity to determine precisely the 
location of the critical point beforehand. 

An estimate of the critical exponent a of the infinite 
three-dimensional Ising system may be obtained by ap- 
plying relation (48) for the maxima of the specific heat of 
the finite cluster systems. Using the system sizes 19 and 
57 one obtains the value a = 0.132, again in reasonable 
agreement with the literature value a = 0.109. 

In Section III we already discussed briefly the possible 
reduction of the correction to scaling terms when going 
from the temperature to the energy picture. Denoting 
the amplitude of the linear correction term in (52) by 

K L = L- a ' v -(B L -b L ) (61) 
a 

the ratio Bl/Kl may serve as a measure of the suppres- 
sion of the linear correction term in (52) compared to the 
corresponding term in (46). Note that the pseudo-critical 
temperature is not integrated into this ratio as T c ^ ap- 
pears in the scaling variable of the energy picture. For 
the magnetisation of the three-dimensional clusters this 
ratio is of the order 8. Similarly the ratio B L /(B L — b L ) 
describes the reduction of the correction term for a sys- 
tem with a logarithmically diverging specific heat. Again 
the factors appearing in the scaling variable of (57) are 
excluded from the ratio. For the two-dimensional Ising 
system in the cluster molecular field approximation this 
ratio is of the order 6. In both cases the corrections to 
scaling are suppressed in the energy picture as expected 
from the considerations in Section III. It is also observed 
in our data that the suppression decreases slightly with 
increasing system size for the range of extensions L in- 
vestigated in our cluster molecular field approximation. 

Finally, we again get a very good estimate of the 
critical temperature from Eq. (37). Inserting the value 
v = 0.63 coming from the best data collapse of the scaled 
data as well as the pseudo-critical temperatures of the 
systems with 7, 19, and 57 spins, we find ksT c /Ji = 
4.5208 which should be compared to the literature value 
/c B T c /J/ = 4.5115. 

C. Three-state Potts models 

Let us close this Section by discussing three-state Potts 
models on square and triangular lattices. In these cases 
the cluster molecular field approximation displays a dis- 
continuous transition for all finite clusters, as illustrated 
in Fig. 3 by the jump of the order parameters. At 



first sight this fact seems to rule out the possible de- 
termination of critical quantities from our small systems. 
However, we show in the following that our methods of 
analysing the data in principle still allow us to obtain 
reliable values for the critical quantities. 
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FIG. 11. Scaled magnetisations of Potts clusters defined on 
a square lattice, see Fig. 3, as function of the reduced tem- 
perature scaling variable. Real data are shown by symbols, 
data obtained from a prolongation of the polynomial fits to 
m sp — 0, see main text, by the full lines. The best data col- 
lapse is achieved for the critical exponent values /3 = 0.11 and 
v = 0.76. 

In the present situation the analysis of effective critical 
exponents has the great advantage that only data com- 
puted for temperatures lower than the critical tempera- 
ture of the infinite system are involved. Therefore the 
discontinuous behaviour of our small systems observed 
at temperatures above the critical temperature does in 
this case not interfere with the determination of the crit- 
ical exponents. Thus the extrapolated critical exponents 
computed for square lattice systems with up to 85 spins 
yield the value (3 — 0.104 for the order parameter crit- 
ical exponent, in good agreement with the exact value 
flex = 1/9. On the other hand the molecular field finite 
size scaling analysis is deeply affected by the disconti- 
nuity of the phase transitions as here data obtained in 
close vicinity to the transition points play a prominent 
role. In fact, the main problem lies in the values of the 
pseudo-critical temperatures which enter the definition 
of the scaling variables. Using the transition tempera- 
tures given by the kinks in the cluster molecular field 
free energies, see inset of Fig. 3, does not lead to a col- 
lapse of the scaled data. A better result is obtained by 
fitting the spontaneous magnetisations of the finite sys- 
tems to polynomials of second order and then inserting 
the temperatures where these polynomial functions hit 
the temperature axis into the scaling variables. This ap- 
proach leads to an acceptable data collapse for the order 
parameter, as shown in Fig. 11, the critical exponents 
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thereby taking the values j3 = 0.11 and v = 0.76. The 
value of (3 again nicely agrees with the exact result but a 
slight discrepancy prevails for v as the exact value is 5/6. 
With regard to the crudeness of our approach, this is still 
a very satisfactory result. However, we did not succeed 
in extracting the critical exponent 7 from the suscepti- 
bility data in a similar way. Finally, inserting the value 
v = 0.76 into Eq. (37) together with the temperatures ob- 
tained from the polynomial fits, we end up with the esti- 
mate ksTc/Jp — 0.998 for the critical temperature, very 
close to the exact value k B T c /J P = l/ln(\/3~+l) w 0.995. 

V. CONCLUSION 

The advantages and shortcomings of molecular field 
theory are well known. On the one hand this approxi- 
mation provides a rather simple way of obtaining ana- 
lytical results even in complicated systems. Usually, the 
results of the molecular field approximation are qualita- 
tively correct. On the other hand, however, large quanti- 
tative discrepancies are found when comparing these pre- 
dictions with exact results or with estimations obtained 
by numerical methods and/or more sophisticated analyt- 
ical approaches. This is especially the case when consid- 
ering systems with a continuous phase transition as the 
molecular field approximation always yields classical val- 
ues for the critical exponents irrespective of, for example, 
the dimensionality of the system. Furthermore, critical 
temperatures are usually overestimated by a considerable 
amount. 

In this work we have shown that the correct values of 
critical quantities can nevertheless be obtained from sys- 
tematic improvements of the molecular field approxima- 
tion. We thereby insert spin clusters of variable sizes into 
a homogeneously magnetised background and determine 
the molecular field acting exclusively on the border spins 
by imposing a condition of homogeneity in the central 
part of the cluster. The interactions inside the clusters 
are treated exactly by computing the density of states 
as function of the internal energy and of the interaction 
energy with the molecular field. 

The data computed in this cluster molecular field ap- 
proximation can then be analysed in various ways in 
order to extract the values of critical quantities. We 
have discussed in this work two rather complementary 
approaches, one approach involving effective critical ex- 
ponents where only data located at temperatures below 
the critical temperature of the infinite system enter, the 
other being a molecular field finite size scaling theory 
(which in the temperature picture is related to that de- 
veloped by Suzuki [6-11]) where the true critical expo- 
nents result from the best collapse of the rescaled data. 
In the latter approach the data located in close vicinity 
to the pseudo-critical temperatures are of main interest. 

With these methods data resulting from the cluster 
molecular field treatment of various two- and three- 



dimensional Ising models as well as of two-dimensional 
three-state Potts models have been investigated. In all 
these cases we have shown that critical temperatures and 
critical exponents can be determined with a surprisingly 
high precision with regard to the smallness of our spin 
clusters. In two dimensions we have the nice situation 
that clusters of various sizes with the same shape can be 
handled by the available computational capacities. This 
is different in three dimensions where our computational 
resources do not allow us to handle spin clusters of vari- 
ous sizes having the same geometry. In fact, these restric- 
tions are imposed by our choice to compute the density 
of states in a numerically exact manner. In principle 
larger system sizes could be investigated by Monte Carlo 
simulations, using one of the highly efficient methods for 
computing the density of states which have been devel- 
oped recently [17]. However, very good numerical data 
would be needed in order to obtain the strength of the 
molecular field as a function of the temperature from 
Eq. (18). Especially the accurate determination of the 
pseudo-critical temperature requires a density of states of 
extremely high precision not easily achievable in a com- 
puter simulation. We therefore expect that data obtained 
from Monte Carlo simulations may not permit a molecu- 
lar field finite size scaling analysis. However, they should 
be very useful for the investigations of effective critical 
exponents, especially when dealing with quantities like 
the susceptibility or the specific heat for which the small 
systems considered in this work did not allow a determi- 
nation of the values of the critical exponents. 

Further progress along the lines developed here can 
be achieved by noting that the homogeneity condition in 
the central part of the cluster renders our cluster molec- 
ular field approximation superior to more traditional ap- 
proaches. However, the spin expectation value is still 
not identical in the whole cluster. Here a scheme devel- 
oped by Galam [12] could be very useful. In that cluster 
molecular field scheme a cluster consisting of the cen- 
tral spin and of its z neighbours is constructed such that 
the homogeneity is extended to the whole cluster. It is, 
however, not immediately obvious how to generalise this 
scheme to larger clusters containing not only the central 
spin and its nearest neighbours. Nevertheless we expect 
that the construction of homogeneously magnetised clus- 
ters (be it by the Galam scheme or by a different one) 
would further improve the systematic cluster molecular 
field approximation and yield even better values for the 
critical quantities. 

APPENDIX 

In this appendix we describe how the density of states 
is evaluated in a numerically exact way. For simplic- 
ity consider a square lattice which consists of L rows of 
length L. In our case the density of states depends on 
two extensive quantities, the interaction energy E of the 
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spins inside the cluster and the interaction energy B of 
border spins with spins a k — 1 outside the cluster. In 
other physical contexts the density of states can depend 
on other extensive quantities, e.g. the interaction energy 
and the magnetisation. The density of states £l(E, B) 
of the cluster is worked out by successively performing 
the trace over the q L states of a row of spins. A simi- 
lar method has already been discussed in the literature 
for a density of states depending only on one extensive 
parameter [24-26] . Let k = 1 , . . . , q L denote the configu- 
rations of the row of spins and I = 1, . . . , L the rows of the 
cluster. Suppose I rows of the cluster are already built 
up and the density of states of the incomplete cluster is 
given by ivi(k, E, B). If the next row is added additional 
internal and external interactions of this row with the 
already built up cluster and the spins outside have to be 
included. Let £k',k denote the interaction energy of the 
microstate k' of the last row of the yet incomplete clus- 
ter with the microstate k of the row to be added in the 
next step. Similarly B k denotes the additional external 
energy. Of course B k does not contain the contributions 
of those bonds of the new row which will become internal 
bonds when the next row is added. Then the density of 
states of the extended cluster consisting now of I + 1 rows 
is given by 

cu l+1 (k,E,B)= ui(k',E',B') 

k',E',B> 

5(E,E> + £ k , jk )S(B,B' + B k ) . (62) 
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The described method is easily generalised to densities of 
states depending on more than two parameters. 
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